Production of spin 3/2 particles from vacuum fluctuations 
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We study the production of spin 3/2 particles in homogeneous scalar and gravitational backgrounds 
using the mode-mixing Bogolyubov method. Considering only the helicity ±3/2 states, we can 
reduce the problem to a standard Dirac fermion calculation and apply the standard techniques in 
a straightforward way. As an example we consider a specific supergravity inflationary model and 
calculate the spectrum of gravitinos created during preheating at the end of inflation. 
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The quantization of fields in the presence of exter- 
nal classical backgrounds leads to interesting phenomena 
such as the production of particles via the amplification 
of vacuum fluctuations. This effect has been mainly stud- 
ied in bosonic models, for example production of scalars 
or gravitons in scalar or gravitational backgrounds. In 
addition, this mechanism for the creation of particles is 
believed to be responsible for the generation of most of 
the particles that constitute the present universe , and 
in fact it plays a key role in the modern theories of pre- 
heating after inflation. In those models, the energy of 
the inflaton field is resonantly converted into particles 
during the period of coherent oscillations after inflation. 
This parametric resonance phenomenon makes the occu- 
pation number of the newly created bosons grow expo- 
nentially fast and causes their spectra to be characterized 
by resonance bands. Recently, the resonant generation of 
spin 1/2 particles has also been considered in the liter- 
ature 1^. In these works, it has been shown that the 
limit on the occupation number imposed by Pauli exclu- 
sion principle is saturated and thus the non-perturbative 
results deviate considerably from what is expected in a 
perturbative approach. 

In this work we are interested in the creation of spin 
3/2 particles through the amplification of vacuum fluctu- 
ations. The generation of such particles in the early uni- 
verse has traditionally been treated by considering the 
perturbative decay of other particles , but not using 
the non-perturbative approach based on the Bogolyubov 
transformations technique. Some estimations of the grav- 
itino production during inflation, based on the analogy 
with Dirac fermions can be found in The spin 1/2 
case suggests that both approaches can give rise to quite 
different results. This could be of the utmost importance 
in the so-called gravitino problem: in supergravity mod- 
els, the superpartner of the graviton field (gravitino) is 
described by a spin 3/2 particle. If such particles are 
created after inflation by some mechanism (particle col- 
lision, vacuum fluctuations) they could disrupt primor- 
dial nucleosynthesis if they do not decay fast enough. 



or if they are stable particles and their masses are high, 
they could overdose the universe. In the perturbative ap- 
proach, these facts impose stringent constraints on both 
the reheating temperature and the gravitino mass . 

The calculation of spin 3/2 particle production from 
vacuum fluctuations is plagued with consistency prob- 
lems that hamper the quantization of such flelds in the 
presence of external backgrounds. It has been known for 
a long time 0, that a spin 3/2 particle in scalar, elec- 
tromagnetic or gravitational backgrounds can give rise, 
apart from algebraic inconsistencies, to faster than light 
propagation modes. This fact completely prevents a con- 
sistent quantization in such cases ^ . The only theory in 
which these problems seem to be absent is supergravity, 
provided the background fields satisfy the corresponding 
equations of motion ||^. However, the complicated form 
of the Rarita-Schwinger equation makes it very difficult 
to extract explicit results even in simple backgrounds. 
In this paper we will show that when we consider helic- 
ity ±3/2 states (which dominate the high-energy inter- 
actions of gravitinos propagating in arbitrary ho- 
mogeneous (and isotropic) scalar or gravitational back- 
grounds, the equations can be reduced to a Dirac- like 
equation. The quantization can be done along the same 
lines as for Dirac spinors and therefore the standard Bo- 
golyubov technique |l^ can be used to calculate the par- 
ticle production. We will also show explicitly, within 
a previously considered supergravity inflationary model, 
that the expected amplification does take place. 

The massive spin 3/2 dynamics in fiat space-time is 
described by the Rarita-Schwinger equation. We will in- 
clude the scalar field coupling by modifying the mass 
term, (following the notation in B): 
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mr,Y]^u^o. (1) 



As usual in supergravity models we will consider Ma- 
jorana spinors satisfying ij]^ — CipJ^ with C = 17^7'^ 
the charge conjugation matrix. Contracting this equa- 
tion with dfj, and 7a7/i we get: 
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(2) 



and 



2^(9A7"V'a- ^V-a) + (^3/2 - $)(7A7''^. + 2Va) = 0. (3) 

Finally contracting this last equation with 7^ we get: 

1(^7'^ V. - 7^ ^iJx) + 3(m3/2 - 'S'W^, = 0. (4) 

When $ = the three equations (||), (||) and (^) can be 
written as the Dirac equation plus two constraints, i.e: 



0, 
0, 
= 0. 



(5) 
(6) 
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The general solution of these equations can be expanded 
in helicity I = s/2 + m modes: 



with Jsm the Clebsch-Gordan coefficients whose values 
are: J-1-1 = Jn — 1, J-n = ^1-1 = l/\/3 and J_io = 
Jio = \f2Ji. Here s) are spinors with definite he- 
licity s = ±f and normalized as u^p,r)u{p,s) — Srs- 
If we set = (Li;,psin0coS(/),psin0sin(/),pcos0) with 
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and p = IpI , then the three spin 1 polariza- 



tion vectors are given by 
1 
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£fj.ip, ^1) = — ^(0,cos6'cos( 
v2 



« sm 0, cos y sm ( 
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9). (11) 



The normalization is e*^(p,m)e^{p,n) — 6mm 
p^ef^{p, m) = p'^e* (p, m) = 0. The corresponding quanti- 
zation details can be found elsewhere Q|. 

Now we turn to the <& ^ case. The expression in (^) 
is no longer a solution of the equations of motion. Let us 
now concentrate on homogeneous scalar fields, only de- 
pendent on the time coordinate $(t). We look for general 
homogeneous solutions of the Rarita-Schwinger equation 
of the form: 



Y] Js,nu{p, s)ef, {p, to) 



s.m 



(12) 



These fields satisfy the condition j^tpfj, — 0, since they 
differ from in just a scalar factor. Now if we re- 
strict ourselves to the helicity I — ±3/2 states, they sat- 
isfy tpQ^^^^ — and, since the spatial derivatives of the 



scalar field vanish {di^ = 0), then (||) and (||) are au- 
tomatically satisfied provided d'^ifii = 0. From (|l^) this 
last condition is equivalent to p^ipi — which holds from 
the condition p^e^(p, m) = 0. Accordingly, for helicity 
±3 /2 states propagating in an homogeneous scalar back- 
ground, the Rarita-Schwinger equation reduces again to 
the Dirac form: 



$(t))v.f/^ 







(13) 



As far as these modes satisfy a Dirac-like equation, it 
appears that all the difficulties in the quantization would 
concern just the helicity ±1/2 modes in this case. In fact 
the above ansatz (|l^) is not a solution for the helicity 
±1/2 modes even for homogeneous backgrounds. 

Let us include the effect of curved space-time. We 
will concentrate on spatially flat Friedmann-Robertson- 
Walker (FRW) metrics, and we will introduce it by 
minimal coupling as done in supergravity, i.e, Upi/'cr = 
{dp + |il^''Eafc)V'cr with Q'^'^ the spin-connection coef- 
ficients and Sah = I [7a, 76]- The e^"''"' removes the 
Christoffel symbols contribution in the covariant deriva- 
tive. We will continue considering ^(t) to be a function 
of time alone. We will only consider the linearized equa- 
tion in 1/M (where Mp — SnM'^) for supergravity [O, 
i.e, we will ignore the torsion contribution to the spin- 
connection which is of 0{M~^). In this case the equa- 
tions of motion for the gravitino read: 

e'^''''"757.^P^^ + ^(^3/2 - <i>)[7^ 7l^. = 0. (14) 

Contracting with Z?^, taking into account that -D^7i/ = 
and [Df^,Dp] = — |i?°^^pSab (the vector part cancels 
because of the e^'^'"^ term), we get: 



- -e^''''"757.^'''^pSab^. - m)Y4'u + (i?^<f )^p 

+ ^(m3/2-<i>)(^7>.-7"W.)-0. (15) 

Following the same steps as in flat space-time we obtain 
(^) and (|^) but replacing ordinary derivatives by covari- 
ant ones. For FRW metrics and helicity ±3/2 states, i.e., 
■00 = 0, it is possible to show that due to the form of the 
Riemann tensor, the first term in (|l5|) is proportional to 
7^-0^ and accordingly we get: 



(tl/)- m3/2 + $)0^ = 0, 



7^V^M 



0, 
= 0. 



(16) 
(17) 
(18) 



Here again we can use the standard formulae for particle 
production obtained for the spin 1/2 cases to study the 
creation of helicity ±3/2 states in a FRW background. 
With that purpose we have to reduce equation (Hq) to a 
second order equation. Let us first write the equation in 
conformal time defined as dt = a{ri)drj: 
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i^ta-'d^:-f-d^ - m3/2 + $ + *2^T°j = 0, (19) 

where a = da/dr]. We will adopt the following ansatz for 
the helicity I = ±3/2 solutions: 



with 



(20) 
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a('7)("^3/2 - *('7))) /p;(?7)u(P, s)'5;iea(p, "i), (21) 



2uj and ^3^2 



and the normalization C^^^(0)C/^(0) 
a(0)m3/2- One can check that this ansatz automati- 
cally satisfies ( |l7|) and ([l8|). An appropriate form for 
the spinor u{j),s) and polarization vectors eaip,n^) can 
be obtained if we choose the Dirac representation for the 
gamma matrices and we take (without loss of general- 
ity) the z-axis to be along the p direction. In this case 
u{p,l)^ = (1,0,0,0), u(p,-l)^ = (0,1,0,0), eaip,l) = 
^(0,l,i,0) and ea(p,-l) = ^(0,1, -i,0). With this 

choice, u{p, ±1) are eigenstates of 7° with eigenvalues 
-|-1. Then equation (|l9[) reduces to the well-known form: 



dif (0(7?) ("i3/2 - *(??))) 

+ a^{vKmy2-Hv)f)fAv)^0- (22) 



In order to quantize the modes we will expand an arbi- 
trary solution with helicity I — ±3/2 as: 



d^p 



(27r)32w 
+ e-^^"^~C/f'^(r/)4-.,) , 



„-3/2(^) (e'P^-[/f(r7)ap, 



(23) 



where the creation and annhilation operators satisfy the 
anticommutation relations {apz,at,p} = (2Tr)^2u}Sii'S{p — 

f)- 

In order to see how this works in practice, we will 
consider a specific supergravity inflationary model (see 
|^,|3|), in which the inflaton field is taken as the scalar 
component of a chiral superfield, and its potential is de- 
rived from the superpotential / = (A^/Af)(0— Af)^. This 
is the simplest choice that satisfies the conditions that su- 
persymmetry remains unbroken in the minimum of the 
potential and that the present cosmological constant is 
zero. The observed CMB anisotropy fixes the inflation- 
ary scale around A = A/Af ~ lO^''. For the sake of 
simplicity, we will consider the case in which the grav- 
itino mass is much smaller than the effective mass of the 
inflaton in this model, 1713/2 ^ ^n^ ~ 10~^M and since 
the production will take place during a few inflaton oscil- 
lations, we will neglect the mass term in the equations. 
The scalar field potential is given by |ll|: 



dl 



'I. 



d(t> Af2 



3Uf 

A-/2 



(24) 



For the above superpotential, the imaginary direction is 
known to be stable and therefore we will take for sim- 
plicity a real inflaton fleld. Along the real direction the 
potential can be written as |Q : 

V{<t>) = A^e*' ((2(0 - 1) + <t>{<t> - Iff - 3(</. ~ \f) (25) 

where we are working in units M = 1. We will assume, 
as indicated in fl^ , that the potential contributions of 
dilaton and moduli flelds are flxed during and after in- 
flation. This potential has a minimum at (/) = 1. The 
coupling of the inflaton field to gravitinos is given by the 
following mass term in the supergravity lagrangian |ll|: 



„G/2 ^ 



If 



(26) 
(27) 



where we have chosen the minimal form for the Kahlcr 
potential G($, <i>^) = <i>^<i> + log j/p. The rest of interac- 
tion terms in the supergravity lagrangian are not relevant 
for our purposes. The inflaton and Friedmann equations 
can be written in conformal time as: 



b ■ 52 
< + 2-(j)+—V 
b 



62 



1 



. 62 
'+A^ 



0, 



V 



(28) 
(29) 



where the derivatives are with respect to the new time 
coordinate 77 = s-^d the new scale factor is de- 

fined as 6(f)) — a{f))/aQ with oq = a(0). The solution of 
this equation shows that after the inflationary phase, the 
scalar field starts oscillating around the minimum of the 
potential with damped amplitude. Substituting in (|2S" 
for this particular case we obtain: 



dr)2 



^^(6e«/^) + ^e")/«K^) = 0, (30) 



with K = p/(aoA'^). From this expression we see that 
when the scalar interaction is switched off, there is no 
particle production, even in the expanding background. 
Following P,E4| we can calculate the occupation number: 



N^i{f) = ^ ( 2k + z[/:,(T)/«,(T) - fimfAT)] 



A2 



6e«(^)/2|A,(T)P 



(31) 



In order for the particle number to be well defined, 
we must evaluate it when the interaction is vanishingly 
small, that is, for large values of T. Here is a solu- 
tion of equation ( ^o|) with initial conditions /^/(O) = 1 
and /k;(0) = —iK which corresponds to a plane wave 



3 



for 77 < 0. In order to define the initial vacuum at 
t) — 0, we have taken the inflaton to be at the mini- 
mum of the potential at that moment (0(0) = 1), which 
implies e^^'^=^'>^^ ~ and b{0) — 1. We have chosen 
4>{0) — 1.8 in our numerical computations which corre- 
sponds to an initial amplitude of the inflaton oscillations 
of about 0.06Mp and a maximum value of the coupling 
e«/2 of 10-i"M„. 
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Figure 1. -Number density of helicity I = ±3/2 
gravitinos (k^A''^;) against k. 

The results for the spectra in the expanding back- 
ground are shown in Fig.l. Note that we have not consid- 
ered the backreaction effects of the produced particles. In 
the flat space calculation, we find that broad resonance 
bands may appear, similar to those in P,p^. When ex- 
pansion is taken into account (Fig.l), the production is 
reduced by 3-4 orders of magnitude, however the number 
of particles produced is not negligible. From Fig. 1, we 
can estimate a lower bound to the total number density 
of gravitinos of both helicities as: 



n{ri) 



1 



TT^a^iv) Jo 



3\6 /-oo 

Npip'^dp = / N^iK^dK (32) 



Tr'^a^ Jq 



Comparing with the number density of a thermal dis- 
tribution of helicity ±1/2 gravitinos as estimated in ||l^] 
(the helicity ±3/2 could be even less abundant) for a typ- 
ical value of the scale factor at the end of inflation ||l^] 
flo ~ 10^^^, the vacuum fluctuation production is sup- 
pressed by a factor 10^°. The corresponding cosmologi- 
cal consequences have been studied in HJl^ . Comparing 
with the entropy density today we obtain: n/s > 10^^^. 
This result is 4-5 orders of magnitude larger than the 
perturbative production during reheating from direct in- 
flaton decay [^3| and it could pose compatibility problems 
with the nucleosynthesis bounds [^J^ for some values of 
the gravitino mass. 

We have considered the production of helicity ±3/2 
gravitinos (which are the relevant states for the cur- 
rent nucleosynthesis bounds) in a particular inflationary 



model. The expression (^2|) shows that the results are 
very sensitive to the model parameters, but they can be 
used to discriminate between different supergravity in- 
flationary models. The completion of the picture would 
require to study other models and also include the pro- 
duction of helicity ±1/2 modes; however, the Bogolyubov 
technique appears very involved for this purpose. (After 
the appearance of this work , the helicity ±1 /2 case was 
considered in Kallosh et al., |iep-th/9907124 ) Acknowl- 
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